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ABSTRACT 

We conducted a survey for infrared excess emission from 16 nearby main sequence shell stars using 
the Multiband Imaging Photometer for Spitzer (MIPS) on the Spitzer Space Telescope. Shell stars are 
early-type stars with narrow absorption lines in their spectra that appear to arise from circumstellar 
(CS) gas. Four of the 16 stars in our survey showed excess emission at 24 /j,m and 70 /xm characteristic 
of cool CS dust and are likely to be edge-on debris disks. Including previously known disks, it appears 
that the fraction of protoplanetary and debris disks among the main sequence shell stars is at least 
48% ± 14%. While dust in debris disks has been extensively studied, relatively little is known about 
their gas content. In the case of (3 Pictoris, extensive observations of gaseous species have provided 
insights into the dynamics of the CS material and surpris es about the composition of the CS gas 
coming from young planetesimals (e.g. IRoberge et al.|[2~006f ). To understand the co-evolution of gas 
and dust through the terrestrial planet formation phase, we need to study the gas in additional debris 
disks. The new debris disk candidates from this Spitzer survey double the number of systems in which 
the gas can be observed right now with sensitive line of sight absorption spectroscopy. 
Subject headings: circumstellar matter - planetary systems: formation - stars: early-type - infrared: 
stars 



1. INTRODUCTION 

Protoplanetary disks may be classified by their rel- 
ative abundance of primordial gas and dust left over 
from star formation. As the amount of this material 
decreases, disks evolve from primordial disks to transi- 
tional disks to debris disks. Debris disks, which are found 
around main sequence stars, contain little or no primor- 
dial material. Instead, they are composed of recently 
produced secondary material generated by collisions be- 
tween and evaporation of planetesimals, analogs of Solar 
System asteroids and comets. Terrestrial planet forma- 
tion may be occurring in the younger debris disks, while 
the older ones seem to correspond to the "clearing out" 
stage early in the history of the Solar System, during 
which most left-over planetesimals are removed from the 
system. Impacts by water-rich planetesimals probably 
delivered most of t he Earth's surface volat iles during the 
clearing phase (e.g. iMorbidelli et al.ll2000l ). 

Gas in debris disks is generally hard to detect. Only 
one bona fide debris disk, 49 Ceti, shows any trace of 
sub- m m CO emission (|Zuckerman et al1 [T995; Den t et all 
2005), indicating that debris disk gas abundances are 
low relative to those in primordial and transitional disks. 
It appears that primordial molecular gas left over from 
star formation has largely dissipated by the debris disk 
phase. However, sensitive UV/optical absorption spec- 
troscopy has revealed small amounts o f secondary gas 
in a few debris disks: Pictoris (e.g. |Lagrange et alj 
119981) . 51 Ophiuchi (e.g. IRoberge et al.ll2002l) . a Herculis 
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(|Chen fc Jurall2003l) . and HD 32297 (|Redfiel J2007I) . The 
gas seen is primarily atomic. The primary production 
mechanisms for secondary gas in debris disks are cur- 
rently unknown, but may include photon-sti mulated des- 
orpti on from circumstellar (CS) du st grains (IChen et al. 
2007]) and/or grain-grain collisions (|Czechowski fc Mannl 
2007). 

(3 Pic is the only debris disk whose gas is well- 
characterized. In this disk, most measured gaseous el- 
ements h^ve_rou£hjy solar composition relative to each 
other (jRoberge et alJ 12006). The exception is carbon, 
which is extremely overabundant relative to every other 
measured element (e.g. C/O = 18 x solar). This discov- 
ery w as surprising, since the central star has solar metal- 
licity (|Holweger et al.lll997h . The carbon overabundance 
suggests that a previously unsuspected - and currently 
unknown - process is operating in the disk, during which 
the planetesimals preferentially lose volatile carbona- 
ceous material. However, the lack of similar abundance 
information for other debris disks renders interpretation 
of the (3 Pic results difficult. 

Detections of gas in debris disks to date have almost 
exclusively used optical or UV absorption spectroscopy 
of edge-on disks. Such observations are sensitive to very 
small amounts of cold gas. Unfortunately, the geometric 
constraint that the line of sight to the central star must 
pass through the disk severely limits the number of disks 
that may be probed for gas in this way. There is a strong 
need for additional debris disks in which gas can be ob- 
served, in order to study the evolution of gas abundance 
and composition throughout the entire planet formation 
phase. 

A class of peculiar stars called shell stars may be the 
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TABLE 1 
Target Stars 





Other Name 


Spectral Type 


TV 

IJistancc 

(P C ) 


Inside trie 
T T"R? 


Exp. Times (sec) 
24 (im 70 (jm 


Main Sequence 


Shell Stars from Hauck & Jaschek (2000) 










HD21620 a 




AOVn 


143 


y 


37 


6965 


HD24863 




A4V 


107 


y 


37 


3105 


HD39283 


£ Aur 


A2V 


71 


y 


37 


671 


HD77190 a 


67 Cnc 


A8Vn 


59 


y 


37 


2685 


HD98353 a 


55 UMa 


A2V 


56 


y 


37 


420 


HD118232 


24 Cvn 


A5V 


58 


y 


37 


294 


HD142926 


4 Her 


B9pV 


148 


y 


37 


2853 


HD158352 




A8V 


63 


y 


37 


797 


HD196724 


29 Vul 


AOV 


65 


y 


37 


671 


HD199603 a 




A9V 


85 


y 


37 


2685 


HD223884 




A5V 


92 


y 


37 


3483 


HD224463 




F2V 


108 


y 


37 


4280 


Shell Stars with 


Time- Variable Inf ailing Gas 












HD42111 a 




A3Vn 


187 


n 


37 


3609 


HD50241 


a Pic 


A7IV 


30 


y 


37 


168 


HD148283 a 


25 Her 


A5V 


79 


y 


37 


671 


HD217782 a 


2 And 


A3Vn 


107 


y 


37 


545 



Binary system. 



avenue toward finding edge-on debris disks containing 
gas. These are stars whose spectra show narrow absorp- 
tion lines arising from line of sight gas at the velocity of 
the star. Most shell stars are B- or A-type. In some cases, 
the narrow lines may arise from interstellar (IS) material, 
although the fact that the gas is at the velocity of the star 
makes this interpretation less likely. In addition, the ab- 
sorption lines are general ly stronger than ca n be readily 
explained by IS gas (e.g. lAbt fc Movdlll973l ). The shell 
stars have high projected rotational velocities compared 
to non-shell early- type stars, showing that gas is more 
likely to be detected when the line of sight intersects 
the r otational plane of the system (e.g. Holwcge r et al.l 
1999). This indicates that the gas is usually confined to 
that plane and strengthens the case for a CS rather than 

15 origin of most shell absorption lines. 

Many shell stars appear to be rapidly rotating evolved 
stars that have recen tly ejected mass from th eir equato- 
rial planes. However. lHauck fe Jaschekl (|2000f ) evaluated 
the luminosity classes of 57 shell stars and found that 
40% of them are main sequence stars. The nature of 
the main sequence shell stars is not fully understood. 
Some are rapidly rotating classical Be stars with hot ex- 
cretion disks. We note, however, that 6 out of t he 23 
main sequence stars in the lHauck fc Jasch ek (2000) sam- 
ple have protoplanetary or debris disks. Among the six 
are the well-known disk systems (3 Pic, HD 163296, and 
MWC 480. In fact, (3 Pic was classified as a shell star a 
decade before discovery of its dust disk fglcttebak 19751 ) . 

In this paper, we present a study aimed at finding addi- 
tional debris disks hiding among the shell stars. We used 
the Multiband Imaging Photometer for Spitzer (MIPS) 
on the Spitzer Space Telescope to measure the fluxes from 

16 nearby main sequence shell stars at 24 and 70 fim. 
The observations were designed to probe sensitively for 
excess flux arising from CS dust. The goal is to assemble 
a unique and valuable set of debris disk systems where 
the gas may be studied in as much detail as the dust. 

2. TARGET SELECTION AND OBSERVATIONS 



We began our target selecti on by considering which 
subset of the 57 shell stars in lHauck fc Jaschekl (2000) 
is the most likely to contain debris disks. To eliminate 
evolved stars, we selected only the stars with luminos- 
ity class V (23 stars). We then eliminated six stars al- 
ready known to have protoplanetary or debris disks. We 
also wanted to eliminate ordinary stars that might have 
been misclassified as shell stars simply because they have 
anomalously large amounts of IS material along their 
lines of sight. Therefore, we discarded stars that are 
outside of the Local Inter stellar Bubble (L I B), us ing the 
LIB size and shape from lLallement et al.l (|2003t ) . This 
left us with 12 stars (listed in Table [T]). 

To this set, we added four shell stars whose 
spectra show narrow, ti me-variable absorption fea- 
tures . These are HD 42111 ( Lccavclicr des Etang s~et al.l 
[1997T). HD 5 241 dHemnel fc SchmittJl2003h. HD 148283 
dGradv et all I1996D. and HD 217782 dCheng fc Neil 
|2003[ ). The time- variability confirms that the gas along 
the lines of sight to these stars is CS, despite the fact 
that one of them is outside the LIB. Such behavior is 
seen in spectra of (3 Pic, and is caused by evaporation 
of star-grazi ng planetesimals p assing through the line 
of sight (e.g. [Beust et al.l fl990). One of the four stars 
(HD 50241) has luminosity class IV. However, since it 
has time- variable infalling CS gas, it is unlikely to be an 
evolved mass-losing star. 

We obtained MIPS 24 /im images and 70 /jm default- 
scale images for 15 of our 16 target stars. For one of our 
stars (HD 77190), we obtained only a 70 ^m observation. 
This star was observed at 24 [im in GTO program 40; 
we used that image in our study. All the data were cali- 
brated with version S13.2.0 of the Spitzer Science Center 
(SSC) calibration pipeline software. The total exposure 
times for each star are listed in the last two columns of 
Table [TJ 

3. ANALYSIS 

3.1. 24 M m Photometry 

To determine the 24 /im fluxes from the target stars, we 
used the pipeline-processed post-Basic Calibrated Data 
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TABLE 2 

Photometry Apertures and Aperture Corrections 



Wave 


'"object 


''sky inner 


'"sky outer 


Aperture 


Comment 


H 


Correction 




24 


14'.' 7 


29'.' 4 


41'.' 7 


1.143 


default a 


24 


14'.' 7 


39" 2 


49'.' 


1.135 


used for HD 42111 b 


24 


14'.' 7 


19'.'6 


29'.' 4 


1.154 


used for HD 158352 b 


24 


9'.' 8 


19'.'6 


31'.'9 


1.322 


used for HD 199603 b 


70 


28'.' 


40'.' 


68'.' 


1.298 


default a 


70 


16'.'0 


40'.' 


64'.' 


1.884 


used for HD 39283, HD 77190, HD 42111, & HD 158352 c 


70 


S'.'O 


40'.' 


64'.' 


4.186 


used for HD 199603 c 



Calculated from our d ata. See ij3.ll 

From the SSC website (http://ssc.spitzer.caltech.edu/mips/apercorr/ ). 



(BCD) mosaic images, which are suitable for detailed 
analysis. The pixels in these images are 2'.'45 x 2'/45. 
Every target star was detected at S/N > 50, except for 
the faintest (HD 224463), which was detected at S/N = 
7. The locations of the star centers were determined by 
fitting a 2-D Gaussian to each star. 

For most of the target stars, we performed aperture 
photometry using an object aperture and sky annulus 
combination close to the 24 /im large aperture combina- 
tion used in ISu et al.1 (|2006l ) (r bj = 14'.'7, r s k y inner = 
29'.'4, r s k y outer = 41'.'7). The aperture correction for 
this combination is 1.143. For three of the targets, 
there are nearby faint stars or background objects that 
fall within the default object aperture or sky annulus. 
In these cases, we used aperture and annulus combina- 
tions that avoid the nearby sources. Aperture correc- 
tions for these combinations were calculated using the 
24 /im image of our brightest target star at this wave- 
length (HD 50241). To verify the accuracy of our cal- 
culated correctio ns, we used tha t image to calculate a 
correction for the ISu et all (|2006f ) aperture combination 
given abo ve. Our c a lculat ed correction differed from the 
published [SueTaD (|2006l ) correction by only 0.5%. All 
the apertures, annuli, and corrections used appear in Ta- 
ble n 

The results of our aperture photometry, which have 
not been color-corrected, appear in Tabled] The uncer- 
tainties in the flux measurements were determined in the 
following way. The aperture photometry equation giving 
the fluxes in mJy is 

F = cA P n£fij-N*-s\ , (1) 

where c is a units conversion factor (0.023504 mJy sr 
arcsec -2 MJy -1 ), A is the aperture correction, p is the 
angular size of the pixels (2.45 2 arcsec 2 ), fij is the sur- 
face brightness value of pixel i,j in units of MJy sr -1 , 
N is the number of pixels in the object aperture, and s 
is the mean sky brightness calculated from the pixels in 
the sky annulus. The sum is over all the pixels in the 
object aperture. This equation may be reduced to 

where C is a constant and B is the estimate of the back- 
ground flux in the object aperture. Propagating uncer- 



tainties gives the following equation for the uncertainty 
in the total flux. 

a F = C* Ij2al >j+ a%, (3) 

where the a Si j are the uncertainties in the surface bright- 
ness values and cs is the uncertainty in the estimate of 
the background flux in the object aperture. 

The <j Si . values were taken from the uncertainty array 
provided by the SSC calibration pipeline. To estimate 
erg, we placed eight apertures and sky annuli in empty 
regions around each star. For each empty aperture, the 
difference between the actual value of B and the value 
estimated from the pixels in the sky annulus was cal- 
culated. The standard deviation of the differences was 
taken as the estimate of as- This method of calculating 
the flux uncertainties includes object noise, noise due to 
variations of the background within the sky annulus, and 
noise due to spatial variations of the background between 
the object aperture and the sky annulus. 

We calculated the median background surface bright- 
ness in each image using pixels between 61" and 135" 
from the target star. The background values range from 
17.3 MJy/sr (low) to 49.7 MJy/sr (medium- high) . The 
definitions of low, medium, and high background levels 
may be found on the SSC website 2 . For our target stars, 
variation of the background between the object aperture 
and the sky annulus appears to be responsible for 18% 
to 67% of the total flux uncertainty. 

3.2. 70 fj,m Photometry 

Since the 70 /im post-BCD mosaic images are not 
suitable for detailed analysis, we used the SSC software 
MOPEX to combine the filtered BCD images. MOPEX 
was set to propagate the uncertainty images from the 
SSC pipeline, and to use multiframe temporal outlier re- 
jection, since all the mosaics have good coverage values 
(> 10 exposures in the regions around the target stars). 
The pixels in the final mosaics are 4"0 x 4"0. One of the 
final 70 ^tm mosaics is shown in Figure [1] Seven of the 
16 target stars were detected at 70 fim with S/N > 3. 
Since many of the stars were not detected, we used the 
telescope pointing information in the image headers to 
determine the locations of the star centers. 

For most of the stars, we used an object aperture 
and sky annulus combination close to the 70 ^m large 

2 http:/ /ssc. spitzer.caltech.edu/obs/bg. html 



4 Roberge & Weinberger 



TABLE 3 
Photometry Results 



ID 




24 fim 




70 fim 






F* -Fmips 
(mjy) (mjy) 


CMIPS ((-^MIPS 
(mjy) 


- F*)/F*) * 100 Excess? F* F M ips o"mips 
(%) (mJy) (mjy) (mjy) 


((Fmips 


-F*)/F*)*100 Excess? 
(%) 



HD21620 


27.31 


34.11 


0.64 


24.88 ±2.34 


yes 


3.3 


24.2 


6.5 


640.1 ± 199.1 


yes 


HD24863 


25.06 


28.01 


0.55 


11.77 ±2.20 


no 


3.0 


< 19.4 


5.3 


< 545.4 


no 


HD39283 


84.22 


84.88 


0.70 


0.79 ± 0.83 


no 


10.1 


< 23.5 


5.7 


< 133.9 


no 


HD77190 


45.08 


48.81 


0.50 


8.29 ± 1.12 


no 


5.4 


< 16.2 


3.5 


< 200.1 


no 


HD98353 


108.95 


114.52 


0.53 


5.11 ±0.49 


no 


13.0 


< 54.1 


12.3 


< 315.3 


no 


HD118232 


111.10 


155.67 


0.49 


40.12 ±0.44 


yes 


13.3 


81.4 


7.6 


513.7 ±57.2 


yes 


HD142926 


28.55 


107.52 


0.51 


276.66 ± 1.79 


yes 


3.4 


25.4 


6.3 


651.6 ± 186.8 


yes 


HD158352 


74.02 


103.44 


0.35 


39.75 ±0.47 


yes 


8.9 


235.1 


7.3 


2551.2 ±82.8 


yes 


HD196724 


73.44 


75.21 


0.46 


2.41 ±0.63 


no 


8.7 


17.2 


4.9 


97.0 ±56.2 


no 


HD199603 


44.71 


50.54 


0.38 


13.04 ±0.85 


no 


5.4 


18.2 


3.8 


239.7 ±71.0 


yes a 


HD223884 


34.45 


37.23 


0.39 


8.06 ± 1.13 


no 


i.l 


< 21.1 


1.5 


< 412.2 


no 


HD224463 


3.79 


4.07 


0.59 


7.50 ± 15.51 


no 


0.5 


< 11.7 


3.6 


< 2469.1 


no 


HD42111 


47.02 


51.13 


0.74 


8.73 ± 1.56 


no 


5.6 


< 9.0 


3.0 


< 59.7 


no 


HD50241 


608.29 


647.56 


0.42 


6.46 ± 0.07 


no 


73.0 


75.9 


9.7 


4.0 ±13.3 


no 


HD148283 


66.95 


71.03 


0.43 


6.10 ±0.64 


no 


8.0 


< 21.2 


5.2 


< 164.4 


no 


HD217782 


88.28 


92.10 


0.56 


4.33 ±0.63 


no 


10.6 


< 41.9 


10.4 


< 296.8 


no 



This "excess" is probably caused by a nearby source. 



aperture combination used in ISu et al.l (|2006[ ) (r Q bj = 
28"0, r sky inner = 40'.'0, r sky outer = 68"0). The de- 
fault object aperture and sky annulus are overlaid on the 
70 /im mosaic shown in Figure [1] For stars with nearby 
faint sources in the 24 /im images, we used aperture and 
annulus combinations that avoided the sources, even if 
they were not visible in the 70 /jm images. The aperture 
corrections for these combinations were taken from tables 
available on the SSC website 3 . All the apertures, an- 
nuli, and corrections used appear in Table O The 70 /im 
photometry results, which have not been color-corrected, 
appear in Tabled 

We calculated the uncertainties in the 70 /im fluxes fol- 
lowing the same procedure used for the 24 /im photome- 
try (described in M3. 1[) . For each star not detected at the 
3<7 level, a conservative upper limit on the 70 /im flux was 
set by adding 3 times the total statistical uncertainty to 
any positive flux measured in the object aperture. The 
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Fig. 1. — MIPS 70 fim image of HD 50241. The mosaic was 
created using MOPEX. The pixels are 4'.'0 X 4'.'0. The default 
object aperture used for the 70 fim photometry is overlaid with 
a solid white line. The default sky annulus is shown with dashed 
white lines. 

3 http:/ /ssc. spitzer.caltech.edu/mips/apercorr/ 



median background surface brightnesses near the target 
stars were estimated from the unfiltered post-BCD mo- 
saic images, using pixels between 48" and 80" from the 
stars. The background values range from 5.8 MJy/sr 
(low) to 16.3 MJy/sr (medium). At 70 /im, background 
variation between the object aperture and the sky annu- 
lus appears to be responsible for 43% to 84% of the total 
flux uncertainty for our target stars. 

4. SPECTRAL ENERGY DISTRIBUTIONS 
4.1. Fitting Photo spheric Models 

To see if any of the target stars show excess IR emis- 
sion, we first needed to determine their expected photo- 
spheric fluxes at 24 /um and 70 /im for comparison to the 
observed fluxes. This was done by fitting stellar mod- 
els to the optical and near-IR fluxes from the stars. B 
and V magnitudes from the Tycho-2 Catalogue (H0g et 
al. 2000) were converted to Johnson magnitudes using 
formulae given in Volume 1 of the Hipparcos and Tycho 
Catalogues 4 , then to flux densities. The near-IR magni- 
tudes in the J, H, and K s bands from the 2MASS Cata- 
log were converted to flux densities using the zero-points 
given in the 2MASS Explanatory Supplement 5 . 

Since all but one of the target stars are within the LIB, 
we did not apply an IS extinction correction to the opti- 
cal and near-IR photometry. Kurucz model photospheric 
spectra, with logg = 4.0 for the stars with spectral types 
earlier than A5 and logg = 4.5 for the cooler stars, were 
fit to the optical and near-IR fluxes using x 2 minimiza- 
tion. The best-fitting model spectrum for each star was 
used to calculate the predicted photospheric fluxes at 
24 /im and 70 /im, which appear in Table [3] 

We then considered whether any of the stars showed 
significant excess flux at either wavelength. The excesses 
were characterized by calculating the percent deviation 
of the measured flux from the predicted photospheric 
flux, while propagating the statistical uncertainties in 
the measured fluxes (see Tabled]). First, we compared 

4 http://www.rssd.esa.int/Hipparcos/CATALOGUE_VOLl/scctl_03.pdf 

5 http: / / www.ipac.caltech.edu / 2mass / releases / allsky/doc/sec6_4a.html 
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the excesses to the uncertainties in the absolute photo- 
metric calibration of MIPS data processed with version 
S13 of the SSC calibration pipeline. The uncertainties 
are 4% for 24 /im data and 7% for 70 fim data (MIPS 
Data Handbook, version 3.2). We conservatively desig- 
nated excesses greater than 5 times these uncertainties 
as significant, which corresponds to > 20% at 24 /jm and 
> 35% at 70 /im. The second criterion for a significant 
excess is that it must be at least 3 times greater than its 
propagated statistical uncertainty. This criterion elimi- 
nated apparent excesses only at 70 /im. 

Four of the 16 stars showed significant 24 /jm excesses 
(HD 21620, HD 118232, HD 142926, and HD 158352). 
These stars also showed significant 70 /im excesses. 
The spectral energy distributions (SEDs) of all the tar- 
get stars are shown in Figure [2 One additional star, 
HD 199603, showed a significant 70 /im excess. However, 
the source visible in the 70 /im image of HD 199603 is 
slightly offset from the expected position of the A-star; it 
is located at the position of a nearby faint source visible 
in the 24 /im image of HD 199603. Therefore, we believe 
that the apparent 70 /im excess from HD 199603 is due to 
contamination from this nearby red source, despite our 
care in aperture selection. The nature of this nearby red 
source is not yet known. 

For the rest of the target stars, the mean 24 /im ex- 
cess is 6.3% ± 3.1%. There is a small positive offset, 
although it is not significant at the 3a level. Ordinar- 
ily, this would indicate a problem with the photosphcric 
modeling or with the aperture photometry itself. How- 
ever, we do not believe this to be case. The photometry 
should be accurate, since the MIPS color corrections for 
early-type stars are extremely small. The photospheric 
modeling appears accurate. The average deviation be- 
tween the fitted effective temperatures and the expected 
effective temperatures (based on the spectral types) is 
only 378 K. The accuracy of the fitted effective temper- 
atures also indicates that IS reddening does not signifi- 
cantly affect the photospheric model fitting. 

The small positive offset measured at 24 /im may be 
real. None of the target stars are unambiguously non- 
excess stars, since they are all main sequence shell stars 
with CS gas. If they are not debris disks or protoplan- 
etary disks, then they may be classical Be stars, which 
show small, ho t mid-IR exces ses from free-free emission 
(see Figure 6 in lSu et ai1l2006D . We plan additional stud- 
ies to clarify the nature of the target shell stars without 
significant MIPS excesses. 

4.2. Fitting Models to Excess Emission 

The next steps were to begin characterization of the 
CS material responsible for the detected excesses and to 
set upper limits on the amount of material around the 
stars without significant excesses. We modeled the total 
SED of each star as the sum of the best-fitting photo- 
sphcric model spectrum and a single-temperature black- 
body dust emission model. The free parameters of the 
total SED model were the blackbody dust temperature 
(Tbb) and the fractional IR luminosity (Lir/T*). The 
model was fit to the MIPS 24 and 70 /im fluxes using \ 2 
minimization. If the star's 12 /im flux from the Infrared 
Astronomical Satellite (IRAS) was a detection, it was 
included in the fitting. IRAS upper limits and detec- 
tions at other wavelengths were not included; the IRAS 



photometry has large upper limits and measurement un- 
certainties compared to the MIPS photometry and did 
not usefully constrain the fit. 

For the four stars with significant IR excess, the dust 
parameters of the best-fitting SED models appear in Ta- 
bled] We estimated the uncertainties in these parameters 
in the following way. We compared total SED models, 
each of which was the sum of a Kurucz photosphcric spec- 
trum model and a single-temperature blackbody emis- 
sion model, to the optical, near-IR, and mid-IR pho- 
tometry simultaneously. All the stellar and blackbody 
model parameters were varied and x 2 values calculated 
over large regular grids of parameter space (which could 
be described as the "brute force" method). This is in 
contrast to the fitting procedure described in the previ- 
ous paragraph, which used the downhill simplex method 
to rapidly find the \ 2 minimum. 

The contours of \ 2 from the brute force minimiza- 
tion were used to determine the Ha uncertainties on 
the model parameters shown in Table 2] This method 
should give accurate estimates of the uncertainties even 
if any model parameter is correlated with any other. The 
contours showed that the stellar parameters (T c g and a 
normalization factor) are correlated with each other, as 
are the blackbody parameters (Tbb and Lm/L+). As 
expected, the stellar parameters are nearly independent 
of the blackbody parameters. This is why the parame- 
ters of the best-fitting model from each brute force \ 2 
minimization agree with those found using the method 
described in the first paragraph of this subsection, in 
which minimization was done first for the stellar model 
and afterwards for the blackbody model. 

We next determined upper limits on Tir/T* for the 
12 stars without significant 24 or 70 /im excesses. The 
procedure was similar to that described in the first para- 
graph of this subsection. For each star, we fit the sum of 
the best-fitting Kurucz model spectrum and blackbody 
emission models with fixed values of Tbb to the IRAS 
12 /im flux if it was not an upper limit and the MIPS 24 
and 70 /im fluxes even if they were upper limits. The % 2 
minimization code used forced the models to 1) stay be- 
low upper limits, 2) stay below Tmips + Ic for detections, 
and 3) stay above the photospheric flux. The best-fitting 
models with fixed values of Tbb between 20 K and 300 K 
were determined. For each of these models, the Tir/L* 
value is the upper limit assuming that particular black- 
body temperature. 

Figure [3] illustrates the general behavior of the Lir/L* 
upper limits as a function of Tbb- Plots of the Tir/T* 
upper limits vs. T BB for HD 217782 and HD 77190 ap- 
pear in the panels on the left. The panels on the right 
show portions of the two stars' SEDs, with the best- 
fitting upper limit SED models for four values of Tbb 
overplotted. The Tir/T* upper limits decrease sharply 
with increasing temperature between 20 K and about 
50 K. This is because we do not have photometry at 
longer wavelengths, which would constrain the fitting of 
SEDs with cold blackbody temperatures. As Tbb in- 
creases, the Tir/T* upper limits increase then decrease, 
forming an absolute maximum (e.g. HD 217782) or local 
maximum (e.g. HD 77190) somewhere between ~ 50 K 
and ~ 150 K. The Tir/T* upper limits then gradually 
increase again. The maximum between ~ 50 K and 
~ 150 K occurs at the temperature of the best-fitting 
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Fig. 2. — SEDs of the target stars. The fluxes in the B and V bands from the Tycho-2 Catalogue and the J, H, and K B bands from 
the 2MASS Catalog are plotted with green diamonds. Color-corrected fluxes from IRAS are plotted with blue circles. The new 24 (im 
and 70 fim MIPS fluxes, which have not been color corrected, are plotted with red squares. Upper limits are indicated with arrows. 
Uncertainties in the flux measurements are overplotted with straight vertical lines. The best-fitting Kurucz model photospheric spectra are 
overplotted with solid black lines. In each panel, the total SED model (the sum of the stellar model and a single-temperature blackbody 
model) is shown with dashed or dotted black lines. The blackbody temperatures and Lir/L* values for these models are given in Table 131 
For the four stars with significant IR excess emission at the MIPS wavelengths (HD 21620, HD 118232, HD 142926, & HD 158352), the 
best-fitting total SED model is shown with a dashed line. For the stars without significant IR excesses, the total SED model showing the 
upper limit on the IR excess emission is overplotted with a dotted line. 
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Fig. 2. — Continued. 
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TABLE 4 
Results from SED Fitting 
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Dust Distance (AU) 
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HD142926 


221.6 ± 19.6 


(8.67 ±0.84) x 10~ 5 


18 


44 


HD158352 


76.1 ±0.79 


(9.29 ± 0.26) x 10~ 5 


62 


191 


HD196724 




< 2.4 x 10~ 6 






HD199603 




< 1.8 x 10~ 5 






HD223884 




< 1.6 x 10~ 5 






HD224463 




< 1.1 x 10~ 4 






HD42111 




< 1.0 x 01~ 5 






HD50241 




< 9.8 x 10~ 6 






HD148283 




< 8.4 x 10~ 6 






HD217782 




< 9.6 x 10" 6 







Note. — The upper limits on Lir/L* are valid for values of Tbb between 50 K and 300 K. 
Calculated assuming blackbody grains. 

Calculated assuming 1 /im-size astronomical silicate grains. 



SED model found if one treats a 70 /zm upper limit as a 
detection (see the right-hand panels of Figure [3]) . 

However, for five of the stars without significant ex- 
cesses (including HD 77190), the largest 2/tr/L* upper 
limit occurs at 300 K (see the bottom panels of Figure[3]). 
In these cases, SED models with fairly high blackbody 
temperatures are consistent with the limits imposed by 
the 12 /im and 24 /im fluxes, indicating that the SED 
fitting is also poorly constrained at shorter mid-IR wave- 
lengths. The results of this analysis are shown in Tabled! 
For the stars without significant IR excess, the Ltr/L* 
upper limit given is valid for all blackbody temperatures 
between 50 K and 300 K. 

5. DISCUSSION 

5.1. Stars With Excesses 

This first survey for debris disks around main sequence 
shell stars was successful; 4 out of the 16 stars we ob- 
served with MIPS show IR excess emission at 24 /im 
and 70 /im. The dust temperatures from simple model- 
ing of their SEDs are cool, similar to those of other de- 
bris disks. These candidate debris disks have best- fitting 
Lia/I* values that are about 30 to 100 times lower than 
that of (3 Pic. The mid-IR excesses of these four stars 
do not resemble those of classical Be stars, whose mid-IR 
fluxes observed with Spitzer show a steep powe r-law de- 
cline with increasing wavelength (see Figure 6 in lSu et al.l 
120061) . 

For the stars with excesses, we calculated the distances 
from the central stars of grains with these temperatures, 
assuming blackbody grains and 1 ^m-size astronomical 
silicate grains. The distance calculated assuming black- 
body grains represents a minimum distance, while the 
one calculated assuming silicate grains is probably more 
accurate. The values appear in Table [4] For HD 21620, 
HD 118232, and HD 158352, the calculated distances 
suggest that their CS dust lies in the Kuiper Belt re- 
gion for an A-star (|Su et al.l l2006f) . The dust around 
HD 142926 is somewhat warmer, and may lie closer to 
the central star. 



We have searched the literature for information on the 
ages of our debris disk candidates. HD 21620 is a mem- 
ber of the a Per open cluster, which indicates that it is 
about 80 Myr old. This age is fairly young compared to 
most debris disks, whose ages range from « 10 Myr (e.g. 
(3 Pic) to more than 1 Gyr. We were not able to find 
any published age estimate for HD 118232. In a future 
paper, we will estimate the age of this star using multiple 
age-dating techniques. 

An estimate of the age of HD 142926 appears in 
IZorec et al.l (|2005). In this paper, the ages of Be stars 
were determined by comparing their positions on the H-R 
diagram to theoretical tracks of the evolution of rapidly 
rotating stars from the zero-age main sequence (ZAMS) 
to the terminal-age main sequence ( TAMS) and beyond 
into the post-main sequence phases. IZorec et al.l (2005) 
estimate the age of HD 142926 to be about 78 Myr old, 
which is w 37% of the star's total main sequence lifetime. 

Based on the IR AS fluxes from HD 158352, 
lOudmaijer et al.l (|1992l ) tentatively suggested that it has 
an IR excess. A new analysis of the IRAS data by 
Moor et al. (2006) reaffirmed this suggestion, which is 
now confirmed by our MIPS photometry. HD 158352 is 
not known to be a member of any young moving group. 
Its age is estimated to be 750^50 Myr, based on the 
star's location on the H-R diagram relative to theoretical 
isochrones (Moor et al. 2006). The IR excess from this 
star is anomalously large given its presumed advanced 
age. In a recent Spitzer MIPS survey for excesses from 
nearby main sequence early-type stars, there are no stars 
older than about 300 Myr that show as large a 24 /im 
excess as HD 158352 (|Su et al.ll200"6l) . 

5.2. Stars Without Excesses 

We here discuss two particular stars from our sur- 
vey that do not have significant IR excesses. HD 50241 
(a Pic) does not show significant excess emission at ei- 
ther MIPS wavelength, although a claim for IR excess 
from t his star has been made based on its fluxes from 
IRAS (|Song et al.ll200lD . This star is often used as a 
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Fig. 3. — Lir/L* upper limits as a function of blackbody temperature. In the left-hand panels, Lir/L* upper limits for HD 217782 and 
HD 77190 are plotted vs. Tbb- The largest Ltr/L* upper limits for 50 K < Tbb < 300 K are indicated with horizontal dashed lines. In 
the right-hand panels, portions of the SEDs of these two stars are shown. The best-fitting upper limit SED models for four values of Tbb 
are overplotted with dotted, dashed, and dot-dashed lines. The color-corrected IRAS 12 fim flux is plo tted with a blue circle. The MIPS 
24 /an and 70 (im fluxes are plotted with red squares; 70 (im upper limits are indicated by arrows. See E|4,2l for a discussion of this figure. 



point-spread function reference for coronagraphic imag - 
ing of the (3 Pic debris disk (e.g. lGolimowski et al.l [2006) . 
so it is fortunate that it does not have a large amount of 
CS dust. 

There has been a claim for star-grazin g planetesimals 
around HD 217782 (ICheng fc Nefl l2003h . which would 
make it a true j3 Pic analog. However, the lack of a sig- 
nificant amount of CS dust near the star argues against 
this scenario. The process that is thought to give rise 
to the large numbers of star-grazing planetesimals in the 
(3 Pic syste m, gravitational pe rturbation by an unseen 
giant planet (jBeust et alJll99C[ ). would also be expected 
to cause dust-producing collisions between planetesimals. 
The time-variable, narrow absorption lines seen in spec- 
tra of HD 217782 seem to require a different explana- 
tion. None of the four stars in our sample with claims 
for v ariable CS gas showe d significant IR excess emis- 
sion. Rcd field et ail IpOOl found the same result using 



Spitzer observations of three other stars with variable CS 
gas. The origin of time- variable, narrow absorption lines 
in spectra of main sequence early-type stars is a topic we 
plan to investigate further. 

5.3. The Disk Fraction Among MS Shell Stars 

We surveyed 12 out of the 23 main sequence shell 
stars compiled in lHauck fc Jaschekl (|2000f ). Six of the 23 
were already known to harbor debris disks or younger 
protoplanetary disks (MWC 480, Pic, HD 97048, 
HD 139614, HD 163296, & HD 179218); these stars were 
observed in other Spitzer programs. The remaining 5 
stars were not observed in our program because they lie 
outside the LIB. However, one of them (HD 144667) is 
an accreting Herbig A e star, which very likely h as a pro- 
toplanetary disk (e.g. iGarcia Lopez et al.ll2006l ). There- 
fore, we consider that there were 7 known prot oplanetary 
or debris disks in the lHauck fc Jaschekl (|2000[ ) sample of 
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main sequence shell stars prior to our Spitzer survey. 

We found four additional disk candidates among the 12 
main sequence shell stars from lHauck fe Jaschekf (2000) 
tha t we observed. T h erefor e, the disk fraction among 
the lHauck fe Jaschekl (|2000l ) main sequence shell stars 
appears to be at least 11 out of 23 (48% ± 14%). We 
have no information about the presence of disks around 
four of the 23 stars; therefore, the disk fraction in this 
set of stars could be as high as 65%. 

At this time, we are not able to properly compare the 
disk fraction among the main sequence shell stars to the 
disk fraction in other sets of early-type stars. The disk 
fraction is higher for young early-type stars than it i s 
for early- type stars with ages > 30 Myr (|Su et al.ll2006f ). 
Therefore, when comparing the disk fraction in different 
sets of stars, the sets should have similar distributions of 
stellar age s. We do not currently kn ow the age distribu- 
tion in the lHauck fc Jasche k ( 2000) main sequence shell 
star set, although we plan to estimate the stellar ages 
using multiple techniques in a future paper. 



6. CONCLUDING REMARKS 

We have found four debris disk candidates in our 
Spitzer survey of 16 main sequence shell stars. These 
stars show IR excess emission characteristic of moder- 
ate amounts of cool CS dust. We have begun follow-up 
observations of all of our shell stars, to more fully char- 
acterize their CS material. 

Based on their SEDs, it seems likely that our candi- 
dates will be confirmed as debris disks. If so, they are 
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likely to be edge-on, which is the most advantageous ori- 
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candidates will likely double the number of debris disks 
with currently detectable CS gas. With this sample, we 
may begin a program to determine why some debris disks 
have both gas and dust and how these different compo- 
nents are related. 

As mentioned in the Introduction, the nature of the 
main sequence shell stars has been something of a puzzle. 
Our Spitzer study suggests that at least half of the stars 
in this class have protoplanetary or debris disks. Surveys 
for shell stars using ground-based optical spectroscopy, 
followed by IR photometry to look for CS dust, may be an 
effective way of finding protoplanetary and debris disks. 
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